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Abstract. Measurements of the polarization of J/ip produced in pion-nucleus collisions are in dis- 
agreement with leading twist QCD prediction where J/ip is observed to have negligible polarization 
whereas theory predicts substantial polarization. We argue that this discrepancy cannot be due to 
poorly known structure functions nor the relative production rates of J/ip and xj- The disagree- 
ment between theory and experiment suggests important higher twist corrections, as has earlier been 
surmised from the anomalous non-factorized nuclear A-dependence of the J/ip cross section. 



INTRODUCTION 

One of the most sensitive tests of the QCD mechanisms for the production of 
heavy quarkonium is the polarization of the J/ip in hadron collisions. In fact, there 
are serious disagreements between leading twist QCD prediction |2|] and experimen- 
tal data on the production cross section of 'direct' J/ip and Xi- We would like 
to advocate that polarization of J/ip provides strong constraints on the production 
mechanisms of J / ip and thus can pinpoint the origin of these disagreements. 

In this paper we will present results on the theoretical calculation of the polar- 
ization of J/ip in ttN collisions. The detailed analysis will be published in a later 
paper [p]]. We found that the polarization of J/ip provides important constraints on 
the nature of the production mechanisms and urge that polarization measurement of 
J/ip should be included in the design of future charm production experiment. 



PRODUCTION RATES OF if) AND X j STATES 

In leading twist QCD, the production of the J/ip at low transverse momentum 
occurs both 'directly' from the gluon fusion subprocess gg — > J/ip + g and indirectly 
via the production and decay of Xi an d Xi states. These states have sizable decay 
branching fractions xi,2 — ► J/ip + 7 of 27% and 13%, respectively. 
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In this model, we assume that the non-perturbative physics, which is described 
by the wave function at the origin in cases of production of J/ip and ip' , is separable 
from the perturbative hard subprocess, i.e., factorization holds. As the wave function 
at the origin can be related to the leptonic decay amplitude, the ratio of ip' to direct 
J /if) production can be expressed in terms of the ratio of their leptonic decay width. 
More precisely, taking into account of the phase space factor, 



a(ip') r(ip'^e+e-) Mj M 

~ ~ 0.24 ±0.03 (1) 



a dir (J/4>) r(J/V e+e") M|, 

where o~dir(J/ip) is the cross section for direct production of the J /if). The ratio ([!]) 
should hold for all beams and targets, independently of the size of the higher twist 
corrections in producing the point-like cc state. The energy should be large enough for 
the bound state to form outside the target. The available data is indeed compatible 
with ([!]). In particular, the E705 values [|J with different projectiles are all consistent 
with 0.24. 

The anomalous nuclear target A-dependence observed for the J/ip is also seen 
for the ip' ||, so that the ratio is indeed independent of A. Therefore, at high 
energies, the quarkonium bound state forms long after the production of the cc pair 
and the formation process is well described by the non-relativistic wavefunction at 
the origin. 

In leading twist and to leading order in a s , J/ip production can be computed from 
the convolution of hard subprocess cross section gg — > J/ipg, gg — > Xj-> e ^ c -> with 
the parton distribution functions in the beam and target. Higher order corrections 
in a s , and relativistic corrections to the charmonium bound states, are unlikely to 
change our qualitative conclusions at moderate Xp. Contributions from direct J/ip 
production, as well as from indirect production via xi an d X2 decays, will be included. 
Due to the small branching fraction xo ~ > J/ip + 7 of 0.7%, the contribution from 
Xo to J/ip production is expected (and observed) to be negligible. Decays from the 
radially excited 2 3 Si state, ip' — > J/ip + X, contribute to the total J/ip rate at the 
few per cent level and will be ignored here. 

In Table 1 we compare the xi production cross section, and the relative rates of 
direct J/ip and Xi production, with the data of E705 and WAll on tt~N collisions at 
Eiab = 300 GeV and 185 GeV Q . The X2 production rate in QCD agrees with the data 
within a 'i^-factor' of order 2 to 3. This is within the theoretical uncertainties arising 
from the J / ip and x wavef unctions, higher order corrections, structure functions, and 
the renormalization scale. A similar factor is found between the lowest-order QCD 
calculation and the data on lepton pair production [jn| On the other hand, Table 
1 shows a considerable discrepancy between the calculated and measured relative 
production rates of direct J/ip and Xii compared to X2 production. A priori we 
would expect the .fT-factors to be roughly similar for all three processes. We conclude 
that leading twist QCD appears to be in conflict with the data on direct J/ip and 
Xi production. Although in Table 1 we have only compared our calculation with the 





a(X2) [nb] 


a dir (J/ip)/a(x2) 


^(xOMxa) 


Experiment 


188 ±30 ±21 


0.54 ±0.11 ±0.10 


0.70 ±0.15 ±0.12 


Theory 


72 


0.19 


0.069 



Table 1: Production cross sections for \ii X2 and directly produced J/ip in n~N 
collisions. The data from Ref. |4], H include measurements at 185 and 300 GeV. The 
theoretical calculation is at 300 GeV. 

E705 and WA11 n~N data, this comparison is representative of the overall situation 
(for a recent comprehensive review see 0). 



POLARIZATION OF THE J /if; 

The polarization of the J/ip is determined by the angular distribution of its decay 
muons in the J/ip rest frame. The angular distribution of massless muons, integrated 
over the azimuthal angle, has the form 

do o „ . . 

oc 1 + A cos 2 6 (2) 



dcos9 



where we take 6 to be the angle between the fi + and the projectile direction (i.e., we 
use the Gottfried-Jackson frame). The parameter A can be calculated from the cc 
production amplitude and the electric dipole approximation of radiative x decays. 

In Fig. la we show the predicted value (solid curve) of the parameter A of Eq. (§) in 
the GJ-frame as a function of xp, separately for the direct J/ip and the xi,2 J/ip+j 
processes. Direct J/ip production gives A ~ 0.25, whereas the production via Xi an d 
X2 result in A ~ —0.15 and 1 respectively. Smearing of the beam parton's transverse 
momentum distribution by a Gaussian function exp [— (fc^/500 MeV) 2 ] (dashed curve) 
has no significant effect in A except for the production via X2 which brings A down 
to A ~ 0.85. The A(x j p)-distribution obtained when both the direct and indirect J/ip 
production processes are taken into account is shown in Fig. lb and is compared with 
the Chicago-Iowa-Princeton J7[ and E537 data || for 252 GeV ttW collisions and 150 
GeV 7i~W collisions respectively. Our QCD calculation gives A ~ 0.5 for xp < 0.6, 
significantly different from the measured value A ~ 0. 

The discrepancies between the calculated and measured values of A is one further 
indication that the standard leading twist processes considered here are not adequate 
for explaining charmonium production. The J/ip polarization is particularly sensi- 
tive to the production mechanisms and allows us to make further conclusions on the 
origin of the disagreements, including the above discrepancies in the relative produc- 
tion cross sections of J/ip, Xi an d X2- If these discrepancies arise from an incorrect 
relative normalization of the various subprocess contributions (e.g., due to higher 



Figure 1: CIP(«) and E537 (o) data compared with theoretical prediction. 



order effects), then we would expect the J/tp polarization to agree with data when 
the relative rates of the subprocesses are adjusted according to the measured cross 
sections of direct J/ip, Xi an d X2 production The lower curve in Fig. lb shows the 
effect of multiplying the partial J/ip cross sections with the required ^-factors. The 
smearing effect is insignificant as shown by the dashed curve. The A parameter is still 
predicted incorrectly over most of the xp range. 

A similar conclusion is reached (within somewhat larger experimental errors) if 
we compare our calculated value for the polarization of direct J / ip production, shown 
in Fig. la, with the measured value of A for ip' production. In analogy to Eq. ([!]), the 
■0' polarization data should agree with the polarization of directly produced J/^'s, 
regardless of the production mechanism. Based on the angular distribution of the 
muons from ip' — > fi + fi~ decays in 253 GeV n~W collisions, Ref. [TT|] quotes X^i = 



0.02 ± 0.14 for xf > 0.25, appreciably smaller than our QCD values for direct J/^'s 
in Fig. la. 



DISCUSSION 

We have seen that the J / ip and Xi hadroproduction cross sections in leading twist 
QCD are at considerable variance with the data, whereas the X2 cross section agrees 
with measurements within a reasonable .ff -factor of 2 to 3. On the other hand, the 
relative rate of ip' and direct J/ip production (Eq. |l|), which at high energies should 
be independent of the production mechanism, is in agreement with experiment. It 
is therefore improbable that the treatment of the cc binding should require large 
corrections. 

In a leading twist description, an incorrect normalization of the charmonium pro- 
duction cross sections can arise from large higher order corrections or uncertainties in 



the parton distributions^. Taking into account that the normalization may be wrong 
by as much as a factor of 10 and that even such a K— factor does not explain the 
polarization data of J/ip, a more likely explanation may be that there are important 
higher- twist contributions to the production of the J/ip and Xi as suggested in large 
xp case |L2], [13 . 

Further theoretical work is needed to establish that the data on direct J / ip and Xi 
production indeed can be described from higher twist mechanisms. Experimentally, 
it is important to check whether the J/ip's produced indirectly via X2 decay are 
transversely polarized. This would show that xi production is dominantly leading 
twist, as we have argued. Thus, the polarization of J/ip production from different 
channels provides a very sensitive discriminant of different production mechanisms. 
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